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Abstract 

The possibility of direct detection of light fermionic dark matter in neutrino detectors is ex- 
plored from a mo del- independent standpoint. We consider all operators of dimension six or lower 
which can contribute to the interaction fp — > e + n, where / is a dark Majorana or Dirac fermion. 
Constraints on these operators are then obtained from the / lifetime and its decays which produce 
visible 7 rays or electrons. We find one operator which would allow fp — ► e + n at interesting rates 
in neutrino detectors, as long as nif < m n . The existing constraints on light dark matter from 
relic density arguments, supernova cooling rates, and big-bang nucleosynthesis are then reviewed. 
We calculate the cross-section for fp — ► e + n in neutrino detectors implied by this operator, and 
find that Super-K can probe the new physics scale A for this interaction up to 0(100 TeV). 



* Electronic address: jenkile@quark.phy.bnl.gov 
< Electronic address: soni@bnl.gov 



I. INTRODUCTION 



Despite years of effort, the identity of dark matter (DM) is still unknown. For a review 
of the evidence for DM, possible DM candidates, and detection techniques, see Here, we 
will be concerned with direct DM detection. Most direct searches for DM consider the case 
where a DM particle scatters elastically off of Standard Model (SM) particles. As DM is 
nonrelativistic, v/c ~ 10~ 3 , the momentum transferred in these interactions is small; for a 
DM particle with mass tudm ~ 100 GeV scattering off a 100 GeV nucleus, the momentum 
transfer is of order 100 MeV. Here, we point out that a similar momentum transfer could 
occur for much lighter DM particles, if they annihilate with SM particles into final-state 
particles with a lower rest mass. Thus, it is reasonable to suspect that rather different mass 
scales could be probed with the same experiments. For illustration, we take the interaction 

fN - f'N' (1) 

where N and N' are SM particles (such as nuclei or nucleons), / is a DM particle (for this 
paper, taken to be fermionic), and /' is some particle which could be either contained in 
the SM or from new physics (NP), but with mf << mj. In this case, the kinetic energy 
of the initial-state / is negligible, but the rest mass mj is converted into kinetic energy of 
the final-state products. If we take rrif to be O(100 MeV) and M N ~ M N > » rrif, the 
momentum transfer can be similar to that of the more standard case given above. 

This thus begs the question of whether or not current experiments could be used to 
rule out DM particles which are very light, but which interact with SM particles via (JT]). 
Recently, there has been much interest in the observability of models with "hidden" sectors 
@, 0, 0,111, where new particles exist with masses similar to those of the SM fields, but whose 
interactions with SM particles are mediated by very heavy particles and, thus, strongly 
suppressed. In this work, we consider a light (here, taken to be below the weak scale) DM 
particle / interacting via the process ([T]) for the case where / is fermionic. We take / to 
belong to a "hidden" sector, i.e., its interactions with SM particles are suppressed by either 
a high mass scale A or by a very small coupling. Although one could consider the case where 
/' is a hidden-sector particle or a neutrino, for simplicity, we consider only dark matter that 
interacts with SM particles through charged- current interactions, i.e., /' is a charged, visible 
SM particle; specifically, for this work, we will assume that /' is an electron. 

Anticipating results which will see in Section HIT] the least-constrained range of DM mass 
is m,f < O(100 MeV). Although we will review the current constraints on MeV-scale DM in 
section 1111} we point out here that DM in this mass range can have an acceptable lifetime 
and relic density and remain compatible with a wide variety of observational constraints. 
The possible range for rrif raises the possibility that the processes fn — > pe~ and fp — > e + n 
could be seen in present-day neutrino detectors which search for neutrinos with energies 
of 10's of MeV. It is the observability of fp — ► e + n which we explore in this paper; we 
will specifically examine constraints obtainable from the search for relic neutrinos from 
supernovae performed by Super-Kamiokande p. 

In this work, we consider the process fp —>■ e + n from a model-independent standpoint 
via effective operators. We do this for two reasons. First, we choose the process fp — > e + n 
based on its experimental observability; as its experimental signature is very similar to SM 
neutrino interactions, this process can be studied with present-day neutrino experiments. 
Second, little is known of the identity (or identities) of dark matter or of possible interactions 
within the "dark sector" or between dark matter and SM particles; these interactions may 
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turn out to be highly complicated and not well- described by any current model. Given these 
considerations, we feel that any plausible method of direct detection of dark matter should 
be investigated, and that a model-independent analysis is in order. 

We point out that our work is somewhat similar in spirit to works on inelastic DM 
@, H, M, EH 11, 12, 13], but with a few important differences. First, the range of DM 



mass which we consider is very different, 0(100 MeV) instead of 0(100 GeV). Second, we 
consider DM which scatters to a lower-mass state; in our case, this lower-mass state is an SM 
particle. 1 Lastly, we point out that in this work, we need not assume that elastic interactions 
such as fp — > fp are suppressed; they could occur, but would have momentum transfers of 
O(100 keV), and thus typical proton recoil energies of O(10 eV), below the energy threshold 
of neutrino experiments. Although the inelastic interaction fp —>■ up could occur with a 
sizeable momentum transfer, the operators which we study in section III 1 1 which contribute 
to this process turn out to be severely constrained. 2 

This paper is arranged as follows. In section [Til we gi ye the basis of operators which can 
contribute to the process fp — > e + n. We derive constraints from DM lifetime in IHH where 
we find that most of the operators are tightly constrained, but that the process fp — > e + n 
can still occur if the DM mass is < O(100 MeV). Given this result, we then review the 
existing constraints on MeV DM. In section Ml we then erive the scattering cross-section 
and give limits on the operator contributing to this process from the Super-K relic supernova 
neutrino search. For the majority of this work, we consider only the case where the SM lepton 
and quark fields are from the first generation. However, in section |Vj we discuss how our 
results change if we consider the same operators, but with different flavors for the quark 
and lepton fields; we also briefly mention issues involved in relating our results to specific 
models. Finally, in section IVIj we conclude. 



II. OPERATOR BASIS 



In this section, we enumerate our operator basis. First, we summarize our notation. We 
denote the lepton and quark doublets of the SM by L and Q, respectively. The corresponding 
SM right-handed singlet fields will be denoted as £r, and d R . We denote the Higgs 
doublet as <f), and define <fr — iT 2 <p*. For now, we take all quark and lepton fields to be from 
the first generation; we delay discussion of other quark and lepton flavors to section [V] 

For simplicity, we consider only fermionic (Majorana or Dirac) DM particles, which we will 
denote as /. We assume that / is a singlet under the SM gauge group, SU (3) x SU (2) x U (l). 3 
We consider operators of dimension six or less which are SU(3) x SU(2) x L^l) -invariant 
and can contribute to the process fN — > f'N' where N, N', and /' are all SM particles. 
Additionally, we require that N, N', and /' are visible SM particles, i.e. not neutrinos. Due 
to the tight constraints on proton decay, we do not consider operators which violate baryon 
number. This limits our current work to only the processes fd — > ue~ and fu — > de + and, 
thus, we concentrate on the process fp — > ne + . Although one could consider processes with 
both visible SM particles and additional invisible, non-SM particles in the final state, we do 



1 DM scattering to a lower-mass state has been considered previously; see 0, 0|; for inelastic, low- mass 
DM in DAMA/Nal, see [16 1. 

2 Here, we do not include operators which contaLrj_bnt 

3 For other studies of fermionic singlet DM, see [l7l . 



j—ppd. non-SM right-handed neutrino fields. 
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not consider that case here. Thus, we consider operators which contain exactly one non-SM 
/ field. We eliminate redundant operators using integration by parts and the equations of 
motion for the fields; although / may have mass or Yukawa-like couplings to other non-SM 
particles, we can still use its equation of motion to eliminate operators containing Jj)j in 
favor of the operators appearing below or operators which do not contribute to the processes 
we are interested in here. 

At dimension four, there is only one operator, L<ftf; this operator contributes at tree level 
to neutrino mass, however, and is thus severely constrained. All operators of dimension five 
can be eliminated by integration by parts and the equations of motion of the fields. Thus, 
we consider operators of dimension six. Again, we ignore operators which are constrained 
at tree level by neutrino mass. With these restrictions, we are left with six operators 



Ow 


= gLr^fW^ 






o VR 


= tmpfuRrfdR, 


o Sd 




Osu 


= Lfu R Q 


O t 





These operators are all suppressed by A 2 , where A is some mass taken to be higher than 
the weak scale, and each operator Oi is multiplied by a coefficient Cj. Note that in all of 
these operators, / is right-handed. Although / has the same quantum numbers as those of 
a right-handed neutrino, we do not, at this stage, require / to have any specific interactions, 
such as a coupling to a right-handed Z 7 , other than those specified above. For simplicity, we 
will neglect the possibility that / interacts through more than one of the above operators 
simultaneously. 

III. DARK MATTER CONSTRAINTS 

Having enumerated the operators Oj which can contribute to fp — > e + n, we now wish 
to place constraints on their coefficients divided by square of the new physics scale, C//A 2 . 
In order for / to be DM, it must be sufficiently long lived; its lifetime must be at least of 
the order of the age of the universe (~ 4 x 10 17 s), and it must decay to visible products 
sufficiently slowly to have not yet been detected. We consider the constraints from DM 
lifetime on each of these operators in turn. 

A. O w 

Ow contributes to the decay / — > 1/7 at tree level; 

tv* \ \C W \ 2 av 2 3 

r(/ -> vi) = — —m f (3) 

where v ~ 246 GeV is the Higgs vacuum expectation value. In order for / to have a lifetime 
comparable to the age of the universe, T < (4 x 10 17 s) _1 = 1.6 x 10~ 42 GeV. However, for 
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/ to be observable via the processes fn — > pe or fp — > e + n at neutrino experiments, its 
mass must be at least a significant fraction of an MeV [2l[ . We obtain 



C w \ 2 . 1 fl MeV 



< 



(4) 

A 4 ~(6x 10 5 TeV) 4 \ m f J K J 

The insistence of a long lifetime for DM is clearly quite constraining. However, an even 

C I 2 I h 

stronger limit on ^4 is possible, using the results of [22], who use limits on 7-ray emission 
from the galactic center determined by INTEGRAL [23[ and measurements of the diffuse 
7-ray background from INTEGRAL Q, COMPTEL [3 and EGRET (SjlSl to constrain 



the lifetime of DM decaying to two daughter particles, one of which is a photon. They find 
that the DM lifetime must be > O(10 26 s) for a DM mass between ~ 1 MeV and ~ 100 GeV, 

with even tighter constraints toward the lower end of this range and down to .04 MeV. Thus, 

\c 1 2 

we obtain a limit on 1 W more than two orders of magnitude stronger than that in Eq. 
dlj), ~ 1/(8 x 10 7 TeV) 4 (l MeV/mf) 3 , a scale clearly inaccessible at neutrino experiments 
or collider experiments such as LHC. 

We note in passing that, although not interesting for current neutrino experiments, masses 



as low as 0.4 keV are allowed for fermionic dark matter [28] . However, limits from [2 



constrain the lifetime of DM particles to be orders of magnitude longer than the age of the 
universe for DM masses as low as 0.2 keV, implying that for even very light O(keV) DM, 
the effects of Ow would be undetectable at both neutrino experiments and at LHC. 



B. Oy 



After Electroweak Symmetry Breaking the Higgs field obtains a vacuum expectation 
value, and O v takes on the effective form 

Cy_ Q ^ -igCyV 2 - 
A 2 v 2V2A 2 

If rrif > 2m e ~ 1 MeV, this allows the decay / — > e + e~u at tree level. For the range of rrif 
valid for the constraints from Super-Kamiokande which we will consider in section IIV} the 
electron and positron masses can be neglected, and O v gives the / a decay width of 

12 1 



ry e+e -„> = ^^.n). ( 6 ) 



Requiring this to be less than 1.6 x 10 42 GeV and taking ntf ~ few MeV, we obtain 

\Cy? . „ / 1 



A 4 ~ ° V(10 3 TeV) 4 J ' ^ 

with even stronger limits for larger values of mf. Thus, we see that the effects of O v are 
out of reach for current-day neutrino and collider experiments for mf > 1 MeV. However, 
even stronger limits are possible; 3(| compute the decay width of an arbitrary DM particle 



to e + e using the 511 keV line measured by INTEGRAL. For the case where / comprises 
all DM, and where its decays are responsible for the 511 keV line, their result corresponds 



to a lifetime t v of 



„ , rt17 10 MeV 
t v ~ 5 x 10 17 yr . (8) 



m f 
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This gives 4 



20 MeV) 

50 MeV) (9) 
80 MeV), 

from which one can see a linear dependence of the NP scale on mj. We note, however, that, 
taken as limits on the new physics scale, these numbers are probably rather conservative. As 
discussed briefly below, [HI and [33[ conclude that the 511 keV line can be produced only 
by processes in which the injection energy of the positrons was less than ~ 3 or ~ 7.5 MeV, 
respectively. Therefore, we note that, for example, if / decays via O v were responsible for 
10% of the 511 keV flux, this would tighten the limits on the scale of new physics A by a 
factor of ~ 1.8. 

We will now use the very tight constraints for O v and Ow to place constraints on the 
remaining operators. 



\C 



— — — ^ ( Tn e 

^4 ~ tn k in5 tw\A4 \""J 



< 



(9.5 x 


10 5 TeV) 4 




1 


(2.4 x 


10 6 TeV) 4 




1 



^ ( 771 f 

~ (3.8 x 10 6 TeV) 4 1 f 



C. O vr 

If rrif > m n + m e , Ovr gives the decay / — > 7i + e~ at tree level. Thus, we consider the 
case where m/ < m^; here, / can decay as / — * e + e~z/ via the mixing of Ovr into O v shown 
in the diagram in Fig. [TJ We note that both quark fields in Ovr are right-handed. Thus, 
the mixing of Ovr into O v must be suppressed by the u and d masses as both quarks must 
flip chirality to interact with the W. The diagram in Fig. [T]is logarithmically divergent; we 
thus estimate the mixing of Ovr into O v as 

Cy{v) _ C VR {A) 1 12m u m d ^A^ ^ 



A 2 A 2 (47r) 2 v 2 m 2 f 



f 



We obtain results of 



\Cvr\ < i 



A 4 ~ (20 TeV) 4 



(m f = 20 MeV) 



< — - (m f = 50 MeV) (11) 

~ (50 TeV) 4 V 1 ' K ' 



< 



' (m f = 80 MeV) 



(80 TeV) 4 



for mf in the range which will be relevant for the limits from Super-K discussed in section 
IIVI Here, like in the case of O v , we see a linear dependence of the NP scale on mf. We 
have taken the values of the light quark masses from [2l[; given the uncertainty in these 
quantities, the error on our estimates of the NP scale should be taken as on the order of 
30%. 



For a slightly weaker bound using decays involving neutrinos, see 31 1. 
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FIG. 1: Mixing of O vr into O v 



However, this expression does not accurately describe contributions from small (< few x 
100 MeV) internal loop momenta. For this region, we instead calculate the contribution from 
a diagram where the internal u and d quarks are replaced by a virtual 7r + . The amplitude 
squared for this diagram is 

spins 

where Gf is the Fermi constant, q = pf — p e , and f w is the pion decay constant, ~ 130 
MeV. Note that this expression contains factors of m/ and m e , analogous to the lepton mass 
dependence that results in the usual SM n + decay helicity suppression. Thus, if the decay 
channel / — > e~yU + z/ M is open, the constraints on CVr/A 2 will be substantially stronger. 
However, for rrif < m^, the contribution to / — > e~e + v from the diagram containing a 
virtual tt + gives the differential cross-section 

dT \C v ^G F \V u tfftmlq\m)-?Y 

dq 2 A 4 10247r 3 m/ {q 2 - ml) 2 [ ' 

which yields the limits 

20 MeV) 

= 50 MeV) (14) 
= 80 MeV). 

These results are somewhat less restrictive than the result obtained via the loop diagram in 
Eq. (1111) . Thus, we take our results in Eq. (TTTi) as our approximte limits on the NP scale. 

Given that Ovr is substantially more weakly constrained than the previous two operators, 
it is fruitful to look for other possible constraints. First, we consider the possible decay 
/ — ► via the mixing of Ovr into Ow We note that Ovr does not mix into Ow at 



VR\ 



A 4 



-G 2 F \V ud \ 2 ftmlm} [ ~ 



q 2 {m 2 — q 2 ) 
(q 2 - ml) 2 



(12) 



\Gvr\ ^ 1 / 

^ [rrif 



A 4 ~ (6TeV) 4 

~ (20 TeV) 4 ijJlf 

^ -, tt {frit 

~ 50 TeV 4 V ; 
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one-loop order. Example diagrams of the contribution to / — > z/7 at two-loop order are 
shown in Fig. [2j however, this contribution is much weaker than the one-loop result for 
/ — > e + e~u, due to an extra suppression from the electron mass; in Fig. [2] (a), this is due 
to an electron chirality flip, while in Fig. [2] (b), it is due to a coupling of the electron to the 
Higgs field. 

We also estimate a limit on the NP scale by the contribution of Ovr to / decay via 
/ — v mixing. This can occur through the diagrams similar to those in Fig. [2J but with 
the external photon removed. The diagram similar to that in Fig. [2] (b) is quadratically 
divergent, and, thus, the contribution to an effective vif mass term is of order 



C VR 2m e m u m d 
(4tt) 4 v 2 



10" 8 eV 



(15) 



as we have cut off the quadratically divergent integral at 0(A) and taken Cvr ~ 1. We 
expect the mixing angle between the / and v to be of order 10~ 8 eV/rrif ~ 1CT 16 . This 
mixing allows the decays / — > vvv and / — > z/e + e~ via couplings to the SM Z or W, 
suppressed by the square of the mixing angle. Thus, we expect the lifetime of the / to be 
of order r M x 10 32 ~ 10 26 s (where is the \i lifetime), longer than the requirement stated 
earlier of ~ 10 17 yr. 

Finally, for comparison, we can also consider constraints from ir decay, as Ovr would 
contribute to 7r + — > e + f. Searches for heavy neutrinos in 7r decay for m/ < 130 MeV 34 
give constraints on Cvr/ A 2 of order 1/(10 TeV) 2 . 





(a) 



(6) 



FIG. 2: Example diagrams for mixing of Ovr, Osd, Os u into Qw at two-loop order. 



D. Sd , Osu, and O t 

Operators Osd and Os u do not mix into Ow or O v at one-loop order. They can be 
inserted into diagrams like that in Fig. Ufa) to give / — > z/7. (Diagrams which contain a 
Higgs, such as in Fig. Wfo), will be suppressed by additional powers of small Yukawa cou- 
plings.) As these operators contain only one right-handed SM field, there are contributions 
with suppression by only one power of a light mass; we give an order-of-magnitude estimate 
for the mixing into Qw'- 

l9 —r- ^ 72 A2 — x 10 • ( 16 ) 



A 2 A 2 (4vr) 4y v \v 2 J A 2 



S 



For rrif = 1 MeV, in order for / to satisfy the 7-ray constraints, we must have 



C w {v) 



< 



1 



(17) 



A 2 ~(8x 10 7 TeV) 2 
with tighter constraints for larger values of rrif. This implies the order-of-magnitude limit 



a 



Su,Sd 



(A) 



A 2 



< O 



(10 3 TeV) 2 



(18) 



Ot-i on the other hand, mixes into Ow at one-loop order and will be more strongly 
constrained than Os u and Osd- We thus find that Os u , Osd, and Ot to be very tightly 
constrained, likely beyond the reach of near-future experiments. 

Given the above constraints, we concentrate for the rest of this work on the operator 
Ovr and restrict ourselves to the mass range rrif < m„ 



It is appropriate at this stage to 
review the motivation for, and the existing constraints on, DM with a mass of order 1 — 100 

its original 
could 
It has 



MeV. Although some supersymmetry- inspired models of MeV DM exist [35j, 136 
motivation has been the proposal [3j 



38|, 1228 that MeV DM annihilating to e+e 
explain the 511-keV line from the galactic bulge observed by INTEGRAL 
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also been argued [42j that searches for 511-keV emissions from DM-rich dwarf spheriodals 



43 1 would provide very compelling evidence for annihilating DM. However, it has been shown 
32| | that in-flight positron annihilations on electrons in the intersteller medium would over- 



produce galactic 7 rays unless their injection energy was < 3 MeV, severely restricting the 
range of DM masses compatible with the 511 keV line. A similar analysis [33| obtained 



an upper limit of 7.5 MeV; see also [4J, |45|, |46|. However, while this makes the case for 



MeV-scale DM less observationally motivated, it merely constrains the annihilation of DM 
to e + e~ at late times, not the existence of MeV-scale DM. Therefore, we will relax the 
requirement that MeV-scale DM be responsible for the observed 511-keV line, and consider 
the entire possible mass range rrif < m n . 

With any DM candidate, its present-day relic density must be compatible with observa- 
tion. The interaction Ovr is similar to that for a right-handed neutrino. As is well-known, 
the Lee- Weinberg bound [47J requires that heavy neutrinos which annihilate through weak- 
scale intermediate gauge bosons have a mass > O(GeV) in order to annhilate sufficiently to 
not overdose the universe. Interactions suppressed by higher physics scales will freezeout 
earlier; if these interactions dominate / annihilations, / will annihilate less completely, thus 
strengthening the Lee- Weinberg bound. Thus, in order for / to be a realistic DM candidate 
with rrif < m n , it must have at least one other interaction in addition to Ovr which is 
stronger than weak interactions. 

A number of works [37J, |38|, [50I . [Bl] . |52| have studied interactions which can make 

the relic density of MeV-scale DM agree with observation. If we assume that / dominates 
the dark matter density, we must have a thermally- averaged cross-section for annihilation 



(cannkrl) at freezeout of O(10~ 25 )cm 3 /s [36j, |37J, |39j, |48j, |50j, |53j, where v r is the relative 
velocity of DM particles. However, if the primary annihilation channel during freezeout 
is // — > e + e~, there is tension between the cross-section implied by the present-day relic 
density and the flux for the 511 keV line, if the cross-section is independent of velocity, i.e., 
s-wave 

E 



36l . l39l . |42| . |48j; similar conclusions hold for gamma rays over a range of energies 
5ll . 54j . However, if the cross-section is p- wave, with (a a n n \v r \) ~ v%, these requirements 
can be brought into better agreement; [38|, 39, 48, 0, Hi, 153 



have suggested obtaining such 
behavior by coupling both DM particles and e + e~ to a new, light U-boson. 1481 find that 
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such a U-boson with an axial-vector coupling to a DM fermion, such as occurs for Majorana 
DM particles, will either give a cross-section that is completely p-wave (for a vector coupling 
to electrons) or whose s-wave component is suppressed by the electron mass (for an axial 
vector coupling to the electron). 5 

Here, as we wish to be model-independent, we will stick with the effective operator 
formalism; for example, for a velocity-dependent annhilation to e + e~, we could introduce 
the operator 

^eVe = (19) 

a a 

where the £r fields are taken to be e flavor. Here, we do not take A a to be similar to the 
NP scale A for Ovr] we take Oy e to be an effective operator valid at energies relevant for 
freezeout, 0(rrif), although one could consider processes mediated by intermediate particles 
which are sufficiently light that this approximation is not valid, if the intermediate particle, 
such as a U-boson, is very weakly coupled. We leave the investigation of how such an operator 
may arise for future work; we make a few brief comments on model-building in section PVT . 
Taking v r ~ 0.4c at freezeout [481 ]. and requiring the thermally-averaged annihilation cross- 
section (er ann |tv|) to be O(10~ 25 cm 3 /s), |CV e |/A 2 is of order l/(few GeV) 2 for / in the mass 
range of interest here 6 . 

However, given the results of 32j and (33| which show that the 511 keV line can only 



be produced by annihilations of DM with very low mass, considering only annihlations to 
e + e~ is not well-motivated. Thus, we also consider the case where / annihilates primarily 
to neutrinos. We can do this by introducing the operator 

^fO Vr = ^h^f R L Tl ,L T (20) 

a a 

or the same operator where the lepton fields are \x flavor. This operator allows the annihi- 
lations // to both t + t~ and v T v T . However, for the mass range of interest here, rrif < 140 
MeV, only the neutrino channel is open both during freezeout and for present-day annihi- 
lations; in the case of the /i-flavored operator, for m/ > m M , annihilations to muons may 



form some fraction of the cross-section [54j. As before, relic density considerations require 
(v a nn\v r \) to be of order 10~ 25 cm 3 /s. Such an operator could again arise from the exchange 
of a light U-boson; however, keeping to our model-independent philosophy, here we again 
only assume it to be some effective operator valid for energies 0(m/). As Oy T does not con- 
tribute to the 511 keV line, a velocity- independent cross-section is acceptable; optionally, we 
could then add in other operators which could contribute to, but not fully explain, the 511 
keV line, as long as they satisfy the requirements imposed by (32I . 33 . 5l| . Upper limits on 
the velocity-averaged annihilation cross-section of DM to neutrinos for DM mass above 65 
and below [66} ~ 100 MeV are few x (10~ 24 cm 3 /s) and O(10~ 25 cm 3 /s), respectively. While 
these limits are not far from the required value for (<T ann |tv|), this annihilation channel is 
not definitively ruled out. Again, taking < cr ann |iv| >~ O(10~ 25 cm 3 /s) gives a physics scale 
for this operator of order several GeV. 



5 There have been many studies of the phenomenolo gy o f a U-boson in high-energy colliders, low-energy 

62 1, neutrino physics [63j |. and its contribution to the electron magnetic 



colliders 15 
moment [5 




For reviews of U-boson detectability, see 



55, 



56|. 



6 We note that if the physics which mediates SM-DM interactions contains very weak couplings, this GeV- 
scale operator could correspond to intermediate particles lighter than a GeV. 
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There are also constraints on MeV-scale DM from supernovae. [67[ placed limits on the 
interactions of right-handed neutrinos by requiring that they interact sufficiently strongly to 
become trapped and not cool the supernova too rapidly, contradicting observations from SN 
1987A. Here, in order to satisfy the relic density requirements, we assume that / interacts 
with electrons and/or neutrinos with a strength much stronger than the weak interactions; 
in this case, the / is trapped. However, if the interaction of the / with neutrinos is too 
strong, it can make the neutrinos more trapped and cause the supernova to cool too slowly; 



681 ] consider the case where a DM particle is strongly coupled to neutrinos and conclude 
that such a DM particle is ruled out if it has a mass < 10 MeV and that this lower bound 
could be raised to as much as ~ 30 MeV. 

The effects of MeV DM on big-bang nucleosynthesis (BBN) have also been studied. 5(3 



consider DM coupled to neutrinos, electrons, and to both neutrinos and electrons simulta- 
neously. They find, for the neutrino case, that DM particles of mass > 10 MeV pose no 
problem for BBN. For the case where DM is coupled to electrons, they find that the cross- 
section necessary for the correct relic density is incompatible with BBN if the interaction 
is s-wave and the DM mass is between 20 and 90 MeV. P-wave interactions, however, were 
still compatible. For the case where the DM particle was coupled to both electrons and 
neutrinos, they expected to have no discrepancy with BBN for a DM mass > 20 MeV. 

Finally, we note that the effects of MeV DM for m/ < 3 MeV on the formation of 
small-scale structure has been studied in j69j . 

Thus, we see that fermionic DM with a mass O(10 — 100 MeV) is not definitively ruled 
out by any of the above considerations, although further investigation would require the 
formulation of a specific model. We now move on to calculating the interaction rate of such 
DM particles in Super-K. 



IV. SIGNALS IN SUPER-K AMIOKANDE 

We now discuss how the / particle could be observed in neutrino experiments. For this 
purpose, we will focus on the relic supernova neutrino search using the process v e p — > ne + 
from Super-Kamiokande 0. We consider this analysis as it covers a wide range of neutrino 
energy E u , and this range in E v can be translated into a range of possible values of m/ in 
the process fp -—>■ ne + . Super-Kamiokande has placed a limit of 1.2 u e /cm 2 s on the relic 
antineutrino flux for 19.3 MeV < E v < 80 MeV. They fit the overall normalization of their 
backgrounds to data, as the energy distributions for their backgrounds and for their signal 
are sufficiently distinct; as the electron in the process fp — > ne + should be approximately 
monoenergetic with an energy of m/ (before being modulated by the detector resolution), we 
assume that this is a conservative estimate of how well this process could be distinguished 
from background. Although it is possible that the difference between the shapes of the 
v e and / energy distributions could result in our obtaining an overly optimistic limit on 
fp — > ne + for specific values of rrif, we would expect this would change the limit on the 
scale of new physics by at most a few tens of percent. 

To calculate the interaction rate of / at Super-K, the first thing which we must know 
is the DM flux. The mass density of DM in our neighborhood is ~ 0.3 GeV/cm 3 0, and 



its velocity Vf relative to the Earth is O(10 3 ) (for c = 1) [711] . Thus, assuming that / 
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constitutes all DM, the / flux <3>j is 



^ 0.3 GeV/cm 3 ^ ^ ^ 10 8 )/cm 2 s (21) 
m/ 

for m/ = (1,10,100) MeV, respectively. Thus, Super-K can be sensitive to the process 
fp — > ne + as long as its cross-section ao is > 10~ 8 times the SM neutrino cross-section 
evaluated at the neutrino energy E u = rrif. The relevant / flux and corresponding sensitivity 
to Co are reduced by a factor of 2 if / and / are present in equal numbers. 

Next, we calculate the ratio of ao to the SM v e cross-section. For simplicity, we neglect 
the electron mass, the proton-neutron mass difference, and the possibility of second-class 
currents. We only include terms leading order in nif/M < 0.1 (and E v /M), where M is the 
nucleon mass. These approximations should not change our values for the NP scale by more 
than O(10%), which is adequate for our purposes. The cross-section ao is 



|2 

167r|v/| A 4 ' 



^ ,.1, J^f mj(\f^ + 3\ 9l \^ (22) 



where /i ~ 1 and g\ ~ —1.27 are the nucelon form factors [72|. For the SM process 



v e p — ► ne + we obtain, using expressions from 72] 



^M^^dAr + SN 2 ). (23) 



7T 



If we take the ratio of these two cross-sections with E v evaluated at mj and insist that this 
be less than 1.2/(cm 2 s)/$j, we find 

|Cv«|V < 1.2/cm 2 s (24) 



8|v/|A 4 ~ (0.3 GeV/cm 3 )|w/|c/m/ 
where, again, v is the Higgs vev and c is in units of cm/s. This gives 

ICvrI 2 „ 1 



< 



A 4 ~ (120 TeV) 



■i 



(m f = 20 MeV) 



< — - (m f = 50 MeV) (25) 

~ (90 TeV) 4 V f ' K ' 



< 



1 (m/ = 80MeV). 



(80 TeV) 



i 



Here, the accessible NP scale is proportional to mj 1 ^ due to the dependence of the DM 
flux on m/. Thus, we see that the reach of neutrino experiments can be quite strong; for 
the mass range considered here, it gives limits on the NP scale of order 100 TeV. 



A few notes about other possibly relevant experiments are in order. First, [73 



uses a 



calculation almost identical to ours to place limits on interactions of right-handed neutrino 



DM using Ga [74|, |75( and CI experiments [76] . However, they considered the mass range 
rrif < 5 MeV, which is dominated by the solar neutrino background; for this mass range, 
the limits on the new physics scale are thus substantially weaker. We also note that limits 
weaker than those above exist on the v P flux for energies below 20 MeV and as low as 1.8 



MeV from Super-K |77j], KamLAND [78j, and the Borexino Counting Test Facility [79]; the 
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limits on the new physics scale for low values of m/ from these experiments would be weaker 
than those obtained above by a factor approximately ranging between a few and an order 
of magnitude. Also, we note that there exist limits on the relic supernova v e (as opposed to 
the v e considered above) flux from the Sudbury Neutrino Observatory j80| and an analysis 
of the Super-K data 81[ which are about two orders of magnitude weaker than that of the 
v e flux. Constraints derived from the v e flux limits are relevant for the cases where 1) / is 
a Dirac fermion and we consider / (and not / as considered above) as our DM candidate, 
and 2) / is Majorana, in which case the flux limits on both u e and v e are constraining. 
We do not consider constraints derived from the v e flux here; however, as the cross-section 
in neutrino detectors is proportional to A -4 , we would naively expect these constraints to 
translate into limits on the NP scale which are weaker than those in Eq. ( [25]) by a factor 
of a few. Finally, we note that facilities such as DUSEL [82| or Hanohano |83j] may improve 
somewhat on these limits in the future. 



V. DISCUSSION 

Having given our main results in section IIV| we now discuss a few issues related to our 
analysis. 

We will first consider how our results change if we consider the same operators, but with 
fields which are not just from the first generation. First, we will consider the effects of 
changing the flavor on the quark fields in operators Ovr, Osd, Os u , and Ot- Since nucleons 
are made of up and down quarks, we must have operators which contain quarks from the 
first generation in order to have contributions to fp — > n£ + . The right-handed quark fields 
in these operators must therefore be of the first generation, while the left-handed quark 
doublet fields can be of the second and third generation, as their lower components contain 
a down quark, suppressed by a Cabibbo-Kobayashi-Maskawa (CKM) matrix term. Thus, 
in the case of Ovr, other flavor structures are not relevant for this work. On the other 
hand, if we change the flavor of the left-handed quark doublet in Os u , the limit which 
we derived by an estimate of the mixing of this operator into 0\y is no longer valid; the 
logarithmically divergent part of Fig. 0(a) is zero by the GIM mechanism. For 0$d, however, 
this diagram is merely CKM-suppressed; the same is true of the one-loop mixing of Ot into 
Ow- It is possible, however, to arrange linear combinations of flavor structures for which 
these constraints no longer apply, but it would be unexpected that the NP scale for the 
unconstrained flavor combinations is substantially lower than that of the flavor diagonal 
case where both quark fields are from the first generation. 

Next, we consider the effects of changing lepton flavor. Our constraints from / lifetime 
for operators Ow, Osd, Os u , and Ot are independent of the lepton mass, except that for p,- 
and r-flavored operators, we are interested in higher values of mf, therefore, changing the 
lepton flavor from e to p or r will only result in stronger constraints. In the case of O v , if 
the lepton is p flavor, / can decay at tree level via / — > p~e + v e as long as m,f > + m e . 
However, as O v is purely charged-current, it can only contribute to neutrino experiments 
where a muon is produced in the final state, which, even neglecting detector thresholds, 
requires mf > m M . A similar argument holds for the r case. Thus, this operator is not 
interesting for neutrino experiments, except perhaps for a very small range of mf, which we 
ignore here. 

Next, we consider Ovr', for simplicity, we will take the lepton field to be p flavor. First, 
we note that / can decay at tree level via / — > p~ir + if rrif > + m n ~ 245 MeV. Thus 
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we are interested in the mass range m M <rrif < 245 MeV. Our expression for the one- loop 
mixing of Ovr into O v in Eq. (TlOT) is independent of lepton flavor, so we expect the /i case to 
be similarly strongly constrained, except for a small mass range rrif ~ m M . However, there 
may be some tension with neutrino mixing, as the mixing term in Eq. ( |15p is dependent 
upon the lepton mass. The observability of / in neutrino experiments for the mass range 
m fj. ~ m f ~ m n + m 7r we leave for future work. We do not address the mass range rrif < 
for this this is not relevant for neutrino experiments. 

Some comments are also in order about the possible applicability of Ovr to specific 
models. Here, we have taken / to be a singlet under the SM. If we assume a single Higgs 
doublet as in the SM, the quantum numbers of the fields in Ovr dictate that / has the 
same quantum numbers as a right-handed neutrino, and that the interaction L(f)f is not 
disallowed by any symmetry. This is problematic, as this operator would allow / to mix 
with neutrinos and decay via / — > vvv. In order to prevent this, we must introduce a second 
Higgs doublet and impose a symmetry which disallows L(f)f. Once this is done, interactions 
must be introduced to give / an appropriate relic density. A complete investigation of these 
points, including whether or not any of the effective operators mentioned in section II III to 
produce a reasonable relic density (or additional operators not considered here, such as 
interactions coupling / to quarks, or interactions involving very light intermediate particles 
where the effective operator formalism is not valid) are generated at an appropriately low 
scale, is beyond the scope of this work. We observe, however, that Ovr appears similar to 
a right-handed neutrino interaction; this may imply that left-right symmetric models could 
be an elightening place to begin constructing a model which contains Ovr as a low-energy 
approximation. 

Finally, another point worth mentioning is that our above limits on the operator coeffi- 
cients assume that / is a DM particle. We have not made any statements on the observ- 
ability of these operators, for example at LHC, if / belongs to a hidden sector, but is not 
a constituent of DM. Thus, the investigation of these interactions, including various flavor 
structures, at LHC could still be a worthy endeavour. 

VI. CONCLUSIONS 

In this work, we have considered the possibility that a light 0(1 — 100 MeV) fermionic 
DM particle / could be observed by direct detection in current-day neutrino experiments if 
it interacts inelastically through the process fp — > ne + . We have approached this problem 
from a model-independent standpoint and constrained operators contributing to this process 
using the DM lifetime and its decays to ju and to e + e~u, finding lower bounds on the new 
physics scale A for most of the operators to be > 1000 TeV. Finding the constraints on 
one operator Ovr substantially weaker than those on the other operators, we then used the 
results from the Super-K relic supernova neutrino search to place limits on Ovr- We find 
that Super-K was able to probe the scale of new physics up to O(100 TeV), far beyond the 
reach of present-day or near-future colliders. 

Given the impressive reach obtainable with Super-K, it is worthwhile to consider possible 
extensions of this work not explored here. First, we have considered only charged-current 
interactions with nuclei; neutral-current interactions such as f\p — > fiV and interactions with 
electrons could also be considered. Also, if neutral-current interactions are considered, it is 
possible that, in addition to neutrino experiments, current direct-detection DM experiments 
designed to look for DM with masses of O(GeV-TeV) could be sensitive to much lighter DM 
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interacting inclastically; the relevance of such interactions to DAMA/Nal has been explored 
in [l6j]. Finally, we have only considered fermionic DM; one could also consider the cases of 
light dark scalars or light dark gauge bosons. 

The main point of our paper, however, is that inelastic interactions in which a DM particle 
annihilates into a particle with a smaller mass (either an SM particle or another as-yet- 
unknown particle) can occur, and that these interactions can be very sensitively probed by 
existing experiments. As we still know very little of the identity, or identities, of DM or how it 
interacts with itself or with "regular" matter, we should consider these possible interactions 
and investigate how existing constraints from neutrino experiments or experiments searching 
for elastically-interacting DM can be translated into constraints on more general scenarios. 
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